The dephasing rate formula in the many body context 
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We suggest a straightforward approach to the calculation of the dephasing rate in a fermionic 
system, which correctly keeps track of the crucial physics of Pauli blocking. Starting from Fermi's 
golden rule, the dephasing rate can be written as an integral over the frequency transferred between 
system and environment, weighted by their respective spectral densities. We show that treating 
the full many-fermion system instead of a single particle automatically enforces the Pauli principle. 
Furthermore, we explain the relation to diagrammatics. Finally, we show how to treat the more 
involved strong-coupling case when interactions appreciably modify the spectra. This is relevant for 
the situation in disordered metals, where screening is important. 



I. INTRODUCTION 

When a quantum system is coupled to an environment, 
one central feature of the resulting dynamics is that the 
quantum system undergoes dephasing, because its de- 
grees of freedom get entangled with those of the envi- 
ronment. Depending on context, a large variety of ap- 
proaches have been developed for calculating the dephas- 
ing rate. In the context of the dephasing of electrons in 
disordered conductors, as measured, e.g., via weak local- 
ization, these include (giving a partial list only): a path 
integral [l[ method to solve the diffusion equation for the 
Cooperon; diagrammatic perturbation theory 0HHH11] ; 
Fermi golden rule (FGR) arguments for the rate of en- 
ergy exchange between the system and the environment 
; an approach relating the loss of phase of prop- 
agating electrons to the change of state of their environ- 
ment dHH or to the loss of purity fl2lfl3l |; a semiclas- 
sical approach [3] using the paradigm of particle plus 
(effective) bath; and elaborations of this idea in terms 
of Feynman-Vernon-type influence functionals for (quan- 
tum)' Nyquist noise . 

These methods vary greatly in their level of rigor 
and/or physical transparency, and in the level of sophisti- 
cation employed in dealing with the subtleties that arise 
due to the indistinguishability of the electron that is be- 
ing dephased from the other electrons that dephase it. 
The associated "Pauli constraints" determine the fate of 
dephasing in the low-temperature limit [3.4.6.7.8,iolfl7l]. 
ensuring that the dephasing rate vanishes in the zero tem- 
perature limit , contrary to some other claims [ToTfltl . 

The present paper provides a pedagogical and physi- 
cally transparent discussion of the role of the Pauli con- 
straints, without undue pretense of generality or rigor. 
We do this within the framework of the so-called SP- 
approach [7.8.9]. Starting from the FGR, it expresses 
the dephasing rate as a dwdq integral over a product 
S(u>, q)P(—tu, —q) that involves two unsymmetrized spec- 
tral functions [IMl9l ]. The first spectral function (S) 
describes the fluctuations of the environment, while the 
other (P) is the power spectrum of the density fluctua- 
tions of the system which characterizes the motion of the 



particles. The clear factorization of the relevant physics 
into system and bath that can exchange frequency (lo) 
and momenta (q) is the main distinguishing feature of 
this approach. 

Previous works using the SP-approach had employed 
a function pW which described the spectrum of a sin- 
gle particle propagating in a fermionic environment. To 
incorporate the physics of Pauli blocking, which plays 
an essential role in determining the upper cutoff on 
the frequency integral, a rather heuristic mix of semi- 
classical and many-body arguments had been employed. 
The present paper aims to rephrase the discussion of 
P in the more general context of an TV-body system. 
The corresponding spectrum can be written down 
using standard and unambiguous many-body construc- 
tions, without recourse to semiclassical arguments, with 
Pauli-blocking factors arising in a very natural and stan- 
dard manner. Remarkably, it turns out that P^] is 
proportional to pW , with the proportionality factor 
Nt = 2T/S F which is the effective number of thermally 
excited particles that can be scattered in a system that 
has a mean level spacing 6 F . This result justifies the 
way in which Pauli-blocking factors had previously been 
built into PW by hand, and places it in a more general 
many-body context. It also clarifies the relation of the 
SP-approach to the influence functional approach [lflfl7T ] 
for dealing with a fermionic system under the influence 
of quantum Nyquist noise. 

The structure of this paper is as follows. Section [TTJ 
briefly reviews the very general "quantum noise ap- 
proach" to the problem of dephasing of a quantum system 
with a discrete spectrum coupled to an environment. Sec- 
tion |TTT] shows that when this approach is generalized to a 
many-body system interacting with an environment, and 
the FGR is invoked to calculate the dephasing rate, one 
readily arrives at an expression of the SP-type, with S 
and P being unsymmetrized spectral functions for system 
and environment. Section [TVl discusses the calculation of 
P for noninteracting fermions in various contexts: P*- -* 
for a single particle; P^ or P^l for an electron or hole 
excitation in a Fermi sea; P' 1 ' for the thermally averaged 
single-particle excitation of the Fermi sea; and P^ for 
the entire many-body system. In Section [V] we use SP- 
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theory to calculate the dephasing rate of a many-body 
system weakly coupled to an environment, and establish 
the remarkable relation = (2T/5 F )P^ mentioned 

above. Finally, Sections I VII and IVlII discuss the case that 
the system and the environment are so strongly coupled 
that screening takes place, which modifies the system 
and/or particle spectral functions, and puts it into per- 
spective relative to diagrammatic approaches. Several 
appendices summarize some technical details, including 
the relation of the SP-approach to a purity-based defini- 
tion of dephasing recently introduced in [T2IIT3I ]. 



belongs to those transitions where the system receives 
some energy. 

It is now a straightforward observation that the total 
decay rate out of some given level can be simplified by 
introducing the spectrum P(u>) of the system operator 
that couples to the fluctuating force. Indeed, we have 



r„ 
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^SMPi-UiK), (4) 



with P(uj;E n ) = FT[(A(t)A{0))] that has the spectral 
decomposition 



II. THE QUANTUM NOISE APPROACH AND 
DEPHASING FOR SIMPLE SYSTEMS 



P(u;;E n )=J2 



A 



2tt5(u - (E m -E n )). (5) 



When considering a quantum system coupled to a dissi- 
pative environment, it is useful to apply the perspective 
of what we term the " quantum noise approach" . This 
means that, at weak coupling, all the effects of the en- 
vironment on the system (dissipation, heating, and de- 
phasing) can be described completely once the frequency 
spectrum of the noisy force coupling to the system is 
known. For a recent general review of the quantum noise 
approach, especially in the context of quantum measure- 
ment and amplification, see [13] • For our specific ap- 
plication to dephasing, we will later employ the scheme 
described in 

In the present section, we will briefly review these ideas 
for the case of a general quantum system with a discrete 
spectrum (e.g. a two-level system). In the following sec- 
tion, we will then extend these considerations, both by 
keeping track of the spatial degree of freedom and by 
preparing the ground for a treatment of many-body sys- 
tems, which is the emphasis of this article. 

Consider a quantum system that couples to some fluc- 
tuating force field JF, such that the interaction is 



V = AT. 



(1) 



The FGR is an expression for the transition rate from an 
eigenstate n to some other eigenstate m of the system, 
and it can be written as: 



r m ^„ = (m An) xS{lu = {E m -E n )), (2) 



Here the quantum noise spectrum S(uo) is defined as the 
Fourier transform (FT) of the autocorrelation function, 



S{u) 



FT 



T(t)T(0) 



T{t)T{0))e wt dt. (3) 



Crucially, such a quantum noise spectrum is asymmet- 
ric, in contrast to the case of a classical stochastic pro- 
cess. This asymmetry contains important physics: The 
weight of the spectrum at positive frequencies indicates 
the rate of processes where energy is released into the 
environment, whereas the weight at negative frequencies 



Note that P obviously depends on the initial system 
state. The structure of Eq. (j4|) nicely indicates that each 
transition corresponds to extracting energy from the bath 
and putting it into the system or vice versa. 

We now turn to the issue of dephasing. First we have 
to agree on a definition for the dephasing rate T v . The 
popular definition is based on a path integral approach 
(see ApplA")). and it has two disadvantages: (i) it be- 
comes ill defined outside of the semiclassical context; 

(ii) it involves an uncontrolled semiclassical (stationary 
phase) approximation which leads to results involving 
symmetrized rather than unsymmetrized spectral func- 
tions, which are thus not consistent with the FGR pic- 
ture. If one does not want to a adopt a context spe- 
cific definition (e.g. relating to magnetoresistance) it is 
advantageous to define T v as the decay rate of the pu- 
rity fl2lfl3T ] (see App |B|) . leading to a result that does 
agree with the heuristic FGR considerations which we 
clarify in the next paragraph, as well as with the more 
sophisticated diagrammatic approach. 

One should be aware that the association between FGR 
transitions and decoherence is not strict for three reasons: 
(i) Different preparations might have different rates of 
decoherence, and consequently there might be (say) two 
rather than one time constants. For example in NMR 
(see e.g.[2l|) there is so called T\ and T% time scales that 
describe the decay of vertical and horizontal components 
of the polarization vector; (ii) A different, non FGR 
mechanism, might be involved. For example in NMR 
the rate 1 /T2 might have a contribution that comes from 
the so called "pure dephasing" type processes which are 
related to energy levels fluctuating in time without induc- 
ing transitions between them. This contribution would 
be given by the fluctuation spectrum at zero frequency; 

(iii) Not any FGR transition implies decoherence, but 
only those that lead to entanglement and hence change 
both the purity of system and that of the environment. 
This is further explained in ApplBl after Eq. (|B2[) . 

With regard to (i) and (ii) we point out that for the 
physical system under study, namely, interacting elec- 
trons in a disordered metal in a thermal preparation, the 
rate is assumed to be well defined in a statistical sense: 
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there is no reason to assume multiple time scales, or the 
existence of a rival mechanism of dephasing that comes 
from zero frequency fluctuations. With regard to (iii) we 
point out that in a more sophisticated treatment, using a 
diagrammatic approach, the elimination of the irrelevant 
transitions is achieved by including "vertex corrections. 
This leads to an effective infrared cut-off in the frequency 
integral ([3]), see for example Eq.(38) or (41) of [lOj. How- 
ever, our main concern here is to elucidate the role of 
Pauli blocking, which will turn out to introduce an ultra- 
violet cutoff into the frequency integral Thus, for the 
purpose of understanding the role of Pauli blocking, it is 
sufficient to ignore vertex corrections and to identify T v 
as T„ of Eq.(j4]), appropriately averaged over the relevant 
energy window as determined by the preparation or the 
temperature. 



III. THE QUANTUM NOISE APPROACH 
APPLIED TO A MANY-BODY SYSTEM 

As we have seen in the previous section, it will be use- 
ful and instructive to write down the formula for the 
dephasing rate of a particle interacting with an environ- 
ment in terms of an integral over a product of spectral 
functions that describe the motion of the particle and the 
fluctuations of the environment. 

In contrast to the preceding discussion, we now want 
to keep track of the spatial degree of freedom explicitly, 
since it becomes relevant in dephasing of particles mov- 
ing in interferometers or a disordered medium. More 
importantly, we also want to consider the general case 
of a many-body system interacting with the environment. 
This will enable us to automatically take into account 
the physics of Pauli blocking which is crucial to correctly 
describe dephasing at low temperatures. The interaction 
between the particle(s) and the environment will be writ- 
ten as 



V 



U{x) p(x) dx . 



(6) 



Here the number density p{x) is either 8{x — x) for a 
single particle, or its many body version in general; while 
U is the fluctuating potential, i.e. an operator which 
is associated with the environment. In the Heisenberg 
(interaction) picture a time index is added so we have 
U(x, t) and p(x, t). 

Following Ref . @Q , we define the spectral functions S 
and P that characterize the fluctuations of the environ- 
ment and the power spectrum of the motion, respectively: 



(ti(x,t)ti(0,0)) 



dtdx, (7) 



(p(x,t)p(0,0)) j**-**" dtdx. (8) 



We assume stationary fluctuations for which the correla- 
tion functions depend only on the time and position dif- 
ferences. Unless otherwise specified the expectation value 



assumes a canonical (thermal) preparation. Note also 
that in Eq. ((SJ) the total volume normalization with L d is 
required in order to get expressions where the infinite vol- 
ume limit is transparent. The spectral function S(q, u>) is 
experimentally well defined: it is essentially the dynamic 
structure factor (note remarks regarding notations in the 
last paragraph of this section). It is measurable in princi- 
ple via scattering experiments, or via the dielectric func- 
tion, to which it is related by the fluctuation-dissipation 
theorem: see Add IU1 and in particular Eq. (|C8p which de- 
scribes the equilibrium fluctuations of the electrostatic 
potential within a dirty metal. Depending on the con- 
text, the physical identification of the spectral function 
P(q, oj) might be more subtle, as discussed in subsequent 
sections. 

We consider a situation where the many-body system 
and the environment are coupled weakly starting at t = 0, 
and calculate the rate T v for transitions induced by the 
coupling. As already explained in the previous section 
we identify this as the dephasing rate. As shown below, 
FGR leads to the following general result [ll|lG3[l3j] , 

/dui ~ 
— S(q,u)P(-q,-u>), (9) 



T v = 



which we call the "SP formula" or the "SP theory" . We 
stress that it is the unsymmetrized (quantum) versions of 
S{q,uj) and P(q,uj) that enter this formula (see 
The purpose of this paper is to discuss the implications 
of this statement in the many body context. 

Equation © can be derived in the standard way. Con- 
sider the probability to have a transition induced by the 
system-environment coupling. To lowest order in the in- 
teraction (i.e. at short times), it reads 



(V(t 2 )V{t 1 ))dt 2 dt 1 



dt%dt 2 / / dx\dx 2 



(10) 



p{x 2 ,t 2 )U(x 2 ,t 2 ) p{x 1 ,ti)U(x 1 ,t 1 ) 

dq—S(q,uj) J J dt\dt 2 J J dx\dx 2 

Employing the standard Golden Rule approach (i.e. go- 
ing to the variables (t\ + t 2 )/2 and r = t 2 — ii, and tak- 
ing the appropriate limit of a total time span much larger 
than the correlation time), we obtain p v ks T v t, with T v 
given above in Eq. ©. 

All the standard caveats of this linear-response treat- 
ment apply. In particular, at very short times, Ptp(t) will 
depend quadratically on time; during intermediate times 
there is a Wigner decay that agrees with FGR; and even- 
tually (latest at the Heisenberg time) the perturbative 
picture breaks down. The condition of weak coupling 
which we impose, is tantamount to demanding that there 
is a large time-window between these two limits, during 
which the FGR- Wigner approximation is valid. 
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It is already apparent from the discussion in the previ- 
ous section that certain conditions have to be met in or- 
der to be able to identify T v of Eq. ([9]) as a dephasing rate 
for some meaningful, experimentally relevant observable. 
The situation we have in mind is that of a particle follow- 
ing different trajectories in an interferometer or traveling 
through a disordered medium. In a diagrammatic lan- 
guage, the loss of interference between these trajectories 
is given by Eq.© provided all induced transitions have 
comparable rates, and provided one is allowed to neglect 
so-called 'vertex corrections' f^liol]. While these vertex 
corrections are indeed important in describing, e.g., de- 
phasing in weak localization, they are not our prime con- 
cern here. 

We end this section with some brief remarks on no- 
tation. For comparison, we indicate the relation to the 
notation adopted in Ref.Q. There it is assumed that 
the potential IA is induced by a background density p. 
The relation between the Fourier components of the po- 
tential and the Fourier components of the background 
density can be written as U q .u = V q ^Pq,u- Assuming 
Coulomb interaction we have V q ^ = Aire 2 /q 2 . Conse- 
quently S(q,u>) = \V qtU \ 2 S s (q,uj), where S s (q,L<j) is iden- 
tified as the dynamic structure factor. For the power 
spectrum of a single particle excitation Rcf. 9] has used 
the analogous notation pW (g, to) = S v (q, lj). These nota- 
tions are oriented for the study of dephasing of electrons 
in dirty metals where the electrons are both the "system" 
and the "bath" at the same time. Using these notations 
the "SP formula" becomes: 



=fdq[ ^\V q ^\ 2 S s (q,Lu)S p (-q,-u J ). 



(12) 



The spectral function S s (q,u>) is experimentally well de- 
fined as explained after Eq.flT]). In contrast, the object 
S p (q,uj) is a theoretical construction, motivated by a 
semiclassical picture (see next section) but not directly 
measurable. In any case both spectral functions repre- 
sent the ability of the bath and the system to exchange 
energy and momentum. Hence the physical meaning of 
Eq. (|12p is quite transparent: It is the sum over all (q,u>) 
modes that allow exchange momentum q and energy lu 
between the particle and the environment. The relative 
minus sign between the (q, u>) arguments of the two spec- 
tral functions reflects the fact that energy (or momen- 
tum) taken by one is given to the other. Thus Eq. (fT2"|) is 
simply the total rate for exchanging " anything" between 
the particle and environment. It can be written down 
almost by inspection. 



PW and P^ N \ describing the dynamics of some ob- 
servable A (e.g. the density) in the context of single 
particle quantum mechanics or many-body Fermi sea 
physics respectively. In particular, we show that the 
many-body spectral function can be written as the 
single-particle spectral function p(°> times the number 
of fermions within the thermally smeared region around 
the Fermi surface. 



A. The single-particle spectrum P^ ' 

We start with a very general discussion of spectral 
functions in non-interacting electronic systems. Given 
any single-particle observable A, we define the single par- 
ticle power spectrum of A in the absence of a Fermi sea 

as 



PW(lo;E) = FT[/i(t)i(0) 



(13) 



^Pn\A mn \ 2 2h8(uj - (E m - E n )), 



Here p n are the microcanonical weights which are peaked 
around E n ~ E, namely p n — S F x S(E n — E), where 5 F 
is the mean level spacing. It is convenient to take the bot- 
tom of the conduction band as the "potential floor" E=0. 
Then if follows from the above spectral decomposition 
that P(co; E) = for to < —E, implying that the poten- 
tial floor provides a lower cutoff for the emission tail. 
If E is well above the potential floor, than the resulting 
spectrum p(°' is essentially classical, i.e. symmetric in to 
provided \u>\ <C E. 

It is implied by the definition Eq.© that the single- 
particle spectral function P(°) (g, lj) is associated with the 
single-particle density operator A = 6(x — x) or equiv- 
alently one may say that the fluctuations of the q 
Fourier component are associated with the special choice 



A = - 



As an important example (to be employed later 



on), we consider the motion of a single particle in a disor- 
dered potential. This motion is diffusive and accordingly 



P<< \q,Lu;E) 



2Dq 2 



(Dq 2 ) 2 



(14) 



where D is the diffusion coefficient. Again we emphasize 
that it is implicitly assumed that the energy E of the 
particle is well above the potential floor, and we regard 
D as a constant within the energy window of interest. 



IV. THE POWER SPECTRUM P(q,v) 

The purpose of the present section is to calculate the 
spectral functions P for the motion of non-interacting 
fermions in an arbitrary system. We distinguish be- 
tween single-particle and many-body spectral functions 



B. The many-body spectrum P l ' 

We now turn to the many-body spectral den- 
sity. If we treat the many body system as a whole 
then we have to employ second quantization to write 
^ = Emn^™n»L«n- Excluding the diagonal n=m 
terms which are irrelevant to FGR transitions we get for 
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FIG. 1: The (q,u>) plane. The power spectra in a metal- 
lic system are distributed pre-dominantly within the shaded 
rectangular area that indicates an implicit momentum cutoff 
(inverse of the mean free path), and an implicit high frequency 
absorption cutoff (related to the rate of collisions). For the 
SP-formula it is essential to realize that the spectral function 
S acl (q, oj) of Eq. (|C8P has a lower emission cutoff which is deter- 
mined by the temperature T. The power spectrum P(—q, — uS) 
of either Eq. (|15fl or Eq. (|21|l . which is associated with the dif- 
fusive motion of a particle, is concentrated pre-dominantly 
within the dark region \uj\ < Dq 2 . For the SP-formula it is 
essential to realize that the energy E of the particle implies 
a frequency cutoff, which is analogous to T. Close to equilib- 
rium one should take E ~ T. 

a non-interacting system in a thermal state 

PM(w) = FT^A(t)A(0)\ (15) 

= ^2 \ A mn\ 2 (aJ.OmO^) 2irS(uj-(E m -E n )) 
nm 

= ^(l-/(£ m ))/(£„)|A mn | 2 27rS(cj-(E m -E n )) 

nm 

— (l-f(E+u))f{E) m(u;E) 




PM(uj;E F ), 



where the last expression is obtained if P<W(w;P) is en- 
ergy independent in the energy range of interest around 
the Fermi energy. Note that, at zero temperature, the 
first factor becomes (5 F T 1 w6'(ijj), cutting off all contribu- 
tions at negative frequencies. Physically, this represents 
the fact that a zero-temperature fcrmion system can only 
absorb energy. If it couples to a zero-temperature envi- 
ronment, there will not be any transitions at all, in con- 
trast to what would be deduced from the single-particle 
spectrum p(°> alone. The overlap between the power 
spectra S and pW is illustrated in Fig[TJ 

It is important to discuss the significance of the many 
body result. To the extent that interactions can be ne- 
glected this result is exact and does not involve any un- 



controlled semiclassical approximation. One realizes that 

P [n] (uj) w -N T xP(°)(w), forw<T, (16) 

where Nt = 2T/S F is the effective number of particles in 
the thermally smeared band around the Fermi surface. 
This factor is extensive, i.e. it grows linearly with the 
system's volume. Only these particles can be excited by 
a small amount of energy u> ~ being absorbed from the 
environment, or, vice versa, they can release some energy 
into the environment. It is crucial to keep the physics 
of Pauli blocking correctly in this description. In fact, 
had we neglected the Pauli blocking, instead of Ea.([l6"|) 
we would have obtained P [N] (u;) = iVP (0) (w), where N 
would have indicated the total number of particles in the 
system. 



C. The single excitation spectrum P l ' 

In the present subsection we would like to make con- 
tact with other descriptions in the literature, where the 
focus is on the dephasing of a single-particle excitation in 
the presence of a Fermi sea. Unlike the many-body cal- 
culation above, such point of view requires to introduce 
Pauli blocking factors " by hand" . One finds (see Eq. [)2ip 
below) that the power spectrum of a thermalized one- 
particle excitation is 

P [1] {lo) « ixP'°)(u), forw<T, (17) 

which implies that Eq. (ll6|) can be re- written as 
pW(u) « N t P [1] (uj), which holds in the vicinity of the 
Fermi energy. One realizes that the factor 1/2 reflects 
that Pauli blocking of the downward transitions and per- 
sists at high temperatures. The relation between the 
spectrum of the single-particle motion and the many- 
body density fluctuations is a central result of the present 
section. It will be employed in the next section to connect 
the one-body and the many-body dephasing rates. 

In the spirit of the prevailing literature we consider 
separately electrons and holes of arbitrary energy, in- 
corporating Pauli blocking factors "by hand" into the 
definition of the particle's power spectrum. For an elec- 
tron above the Fermi sea it has been claimed in Ref . @@| 
that the Fermi energy Ef is like an effective potential 
floor. This is implied by the Fermi statistics, taking into 
account that a one-body-operator can change the state 
of only one electron. It also results from a more de- 
tailed analysis of the consequences of Pauli blocking on 
dephasing, which has been carried out both in the con- 
text of weak localization [6lfic| and ballistic interferom- 
eters [22j|23|, and which has been illustrated in exactly 
solvable models as well [24H251 ]. Consequently we define 
the power spectrum of an electron excitation with the 
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appropriate Pauli blocking factors incorporated: 
P [e \uj-E) = ^2(l-f(E m ))p n \A mn \ 2 2TrS(uJ-(E m -E n )) 

urn 

= (l-f(E+co)) PW(u;E). (18) 
An analogous expression is introduced for holes: 
P [h \w-E) = J2f( E m)Pn\A mn \ 2 2n5{LJ~(E n -E m )) 

urn 

= f(E-cu) P(°'( W ;£). (19) 
We can thermally average over E using the prescription 

G{E) [-f{E)]dE 



G(E) 



(20) 



where f(E) is the Fermi occupation function which is de- 
termined by the Fermi energy E F and the temperature T. 
Then we get 



PW(w) = P^{lo-E) = (21) 
d 



du> 



1 



o/T 



xPW(uj;E F 



Note that, at zero temperature, the first factor is just 
a step-function 0(w), cutting off the contributions from 
negative frequencies. 



D. Density fluctuations 

By specializing the above general discussion to the 
case that A represents the density operator, and us- 
ing the reasoning of the present section, we can define 
the Pi 6 ' (q, uj; E) of an electron in the Fermi sea; the 
P^(q, uj; E) of a hole in the Fermi sea; the P^(q, to) 
of a single particle excitation at equilibrium; and the 
pW(g,w) of the whole many body electronic system. 

The results Eq.®, Eq.(TH]), Eq.® and Eq.®, as 
well as the approximations Eq. (fT7|) and Eq. (fT6|) for the 
power spectrum of the system are all of the form 



P(q,u) = f p (u)P^(q,cj). 



(22) 



where f p (uj) reflects the way in which the Fermi occupa- 
tion statistics manifests itself. 



V. DEPHASING OF A MANY-FERMION 
SYSTEM FOR WEAK COUPLING TO THE 
ENVIRONMENT 

The purpose of this section is to work out the many- 
body dephasing rate r[f ' according to SP-theory, and 
to compare it with the single-particle dephasing rate 
that is obtained after incorporating Pauli blocking factors 
into the spectra, as described in the previous section. 



We will see they are equal, up to a factor describing the 
number of thermally excited particles that effectively can 
participate in the processes leading to dephasing. 

Before we derive the dephasing rate, we set up a few 
simplifications in the notation. The integrand of the "SP 
formula" Eq. ^ includes a product of two spectral func- 
tions. We assume that the "bath" is in thermal equilib- 
rium and therefore the detailed balance condition that 
should be satisfies is: 



o/T 



S(q, -uj) 

It follows that the spectral function can be written as 

2uj 



(23) 



S(q,u>) 



1 - e^/ T 



v(q,u), 



(24) 



where rj(q,uj) is a symmetric function. It represents a 
generalized friction coefficient (in analogy to the standard 
notation in the Caldeira-Leggett model). Note that at 
high temperatures we have S(q,u>) = 2rjT. 

Without any approximation involved, because the du> 
integration extends from — oo to oo, the integrand of the 
"SP formula" can be symmetrized using the replacement 

F(uj) -> {F( w )}^=i[F(w) + F(-u;)]. (25) 



It is now natural to combine the prefactors in Eq. 
and Eq. (|2"2"|) into a "frequency-dependent temperature" 
for u> transitions: 



2u> 



1 - c-"/ T 



(26) 



and to define an associated symmetrized spectral func- 
tion for the effective thermal fluctuation of the environ- 
ment: 



5(0) ( 



q,u) 



2 V (q,io)T(io). 



(27) 



These definitions allow the "SP formula" to be written in 
a symmetrized form that involves a product of functions 
that are symmetric in uj. Namely, 



T v = \JdqJ ^C)( g , W )p(»)( 9|W ). 



(28) 



In the following subsections we shall discuss the func- 
tional form of T(u>) which is crucial for the calculation 
of low temperature dephasing. But first let us illumi- 
nate the outcome of Eq. (f2"51) in what can be termed the 
semiclassical Nyquist limit. Namely, considering high 
temperatures, for which not only p( ' but also are 
classical-alike, one realizes that Eq. (f2"5)l still contains a 
non-classical 1/2 due to the Pauli blocking of the down- 
ward transitions. So strictly speaking Eq. ([28|) does not 
possess a classical limit. If t] is independent of ui and 
the motion of the particle involves only small frequencies 
\u>\ < T, then T^(uj) in the integrand of Eq.([2"8"]l can be 
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replaced by the temperature T, and we get the simple 
result 



rW = aT, 



(29) 



where the dimensionless a is the dq integral over 77. But 
if we consider (say) a diffusive electron, then at low 
temperatures its power spectrum is broader than T if 
q > {T/D) 1 / 2 , as illustrated in Figffl Then the weight 
which is provided by T(lo) is like an effective cutoff, lead- 
ing to non-linear dependence on the temperature. See 
e.g. 0Q for a gallery of various results. 



A. The dephasing rate of quasi particles at 
equilibrium 

The calculation of the one body dephasing rate T$ 
involves the spectral function P^(q, u) of Eq. (|2"Tj) and 
hence 



{ 


2lo 


d 


to 


} 






dtu 


1 -e"/ T _ 





(30) 



The calculation of the N body dephasing rate in- hence 
volves the spectral function P^ N '(q,u>) of Eq. (fT5|) and 
hence 



T( W ) 



[JV] 



2w 



—u>/5 F 



Doing the algebra we get 

(w/2T) 



SI/ ill 



T(uj 



1 2 



_sinh(w/2T) 
and T(u)W = (2T/<5 F )T(w)W, leading to 



r. 



(31) 



(32) 



(33) 



We have thus come to the conclusion that the many- 
body dephasing rate is equal to the single-excitation 
dephasing rate T$ which properly incorporates Pauli 
blocking, multiplied by a factor Nt = 2T/S F that counts 
the number of thermally excited particles. At this point 
it is important to emphasize that if the interfering entity 
(system) consists of N s constituent particles, all interact- 
ing in the same way with the environment, one expects 
the dephasing rate to be N s times the dephasing rate 
for a single constituent particle under the same condi- 
tions. This is why larger systems are expected to deco- 
here faster (26||27| . In our case the effective number of 
participating particles in the system is Nt irrespective 
of the actual total number. Thus, without putting in 
"by hand" any Pauli blocking factors, we have re-derived 
the correct result that is obtained from a more sophis- 
ticated diagrammatic or path-integral analysis [6tfl0j|. It 

is therefore possible to regard v\p as the dephasing rate 
per effective particle. Equation (|3"3"|) represents a cen- 
tral result of the present paper, whose consequences and 
modifications will be discussed in the following. 



B. The dephasing rate of a non-equilibrium 
excitation 

It is interesting as well to consider the dephasing rate 
for a non-equilibrium one-particle excitation at some en- 
ergy E. Then we cannot simply start from the many- 
body spectrum P^ N ' which is calculated for the thermal 
equilibrium state of the many-body problem. To make 
contact with previous approaches in the literature [ic| |. 
we briefly formulate the nonequilibrium dephasing rate 
in terms of the present notations. In Ref.[10( it has been 
argued that the dephasing rate should be calculated as 



^E)^ l -{T^{E)+Y^{E) 



(34) 



If we make a thermal average over E both terms are 
equal, but if we consider non-equilibrium excitations a 
more careful treatment is required. Using the expres- 
sions for pl e l(q,w; E) and P^l (q, lo; E), the correspond- 
ing function in Eq. (f2"2")) is 



f p (tu) - [1 - f(E+to) + f(E-u;)}/2, 



(35) 



T(u,) 



l_ e -u/T 

Doing the algebra we get 
T(w) 



(l-f(E-uj)+f(E+ujj) f (3(i! 

sym 



( E+u Y 



, / U3 \ 1 , / E—LU 

coth — H — tanh 

\2T J 2 \ 2T 



tanh 



u>, 



which agrees with Ref . [lj| • 



VI. DEPHASING OF ELECTRONS WITH 
SCREENING: THE DIAGRAMMATIC 
PERSPECTIVE 

Up to now, we have treated the weak-coupling case, 
in which neither the particle dynamics (i.e. the density 
fluctuations) nor the potential fluctuations are influenced 
appreciably by the presence of interactions. It is for this 
reason that we have been able to derive the total particle 
decay rate in such a simple manner, expressing it through 
the product of correlators P and S, associated with the 
particle dynamics and the potential fluctuations, respec- 
tively. 

More care has to be exercised when one tries to extend 
the framework to cases in which the potential fluctua- 
tions themselves change strongly (even qualitatively) af- 
ter switching on the interaction. The most important ex- 
ample concerns the effect of the long-range Coulomb in- 
teraction in a metal. The long-range repulsion suppresses 
efficiently long-wavelength, low-frequency density fluctu- 
ations and the associated potential fluctuations, thus giv- 
ing rise to screening. An important counter-intuitive (yet 
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(a) 



(c) 



(b) 



N' 



(d) 




N 



(e) 



1 




FIG. 2: The relation between "SP theory" and standard 
diagrammatics: (a) Diagrammatic calculation of the single- 
particle decay rate to leading order in the density-density in- 
teraction, (b) When applying the "SP theory" to a many- 
body system, one effectively calculates an extensive diagram, 
whose value scales with the number of affected particles (here 
denoted as "N"). (c) For a metal with screening, diagrams 
of this type have to be summed to yield the single-particle 
decay rate in RPA. (d) After translation into the diagrams 
appearing in the "SP theory", we see that the bare particle- 
hole bubble has to be used in describing the system's motion, 
in contrast to (e). 



well-known) consequence is that the potential spectrum 
at small ui,q is independent of the electron charge (see 
Eq.([CH|). 

It will turn out that, for such a situation, the single- 
particle decay rate is not correctly reproduced with a 
naive ansatz, in which both S and P are obtained for the 
full interacting system (neither, of course, can we neglect 
interactions completely in both S and P). To shed light 
on this issue, we now rephrase the results of the "SP 
theory" developed here in terms of diagrams. 



A. Diagrams for the weak-coupling limit 

We first return to the weak-coupling case. In that 
limit, the single-particle decay rate is obtained in a di- 
agram of the type shown in Fig[2]Ja). It represents 
the interaction of the given particle with the density- 
fluctuations described by the bubble. Turning to SP- 
theory, we see the following: Its straightforward applica- 
tion to a many-body system produces a diagram of the 
type shown in FigfS^b). This diagram has no external 
lines, and is therefore extensive, i.e. its value grows with 
volume. More precisely, as elaborated in our previous dis- 
cussion, it yields the decay rate multiplied by the number 
of particles that can be scattered. Once this feature is 
taken into account, one can deduce the single-particle 
decay rate. 



electron(s) 


bath of phonons 




|electron(s) 


the other 
electrons 




Selections 

FIG. 3: We first consider non-interacting electrons in a box 28 
of volume L d coupled to a bath of phonons at temperature T. 
The solid curve on the right illustrates how the decay rate of 
the purity depends on the number of electrons in the system. 
If the electrons were not identical (dotted line) one would get 
r = AT' 1 ', but due to the Fermi statistics one obtains a sat- 
uration at r = iVrr' 1 '. In the case of N interacting electrons 
in a dirty metal, with no bath of phonons involved, one can 
define a bunch of N s test particles as the system^, while all 
the other N—N 3 electrons constitute the fluctuating environ- 
ment. We have in mind N s < N, while for (N-N s ) < N, 
the role of system and environment is flipped. 



B. Diagrams beyond weak-coupling 

Let us now have a look at the strong-coupling case, 
where screening alters drastically the fluctuations of the 
density and the potential. Within the random phase ap- 
proximation (RPA) , the single-particle decay rate can be 
calculated using diagrams of the type shown in FigJ^c), 
with an arbitrary number of polarization bubbles inserted 
to account for screening. In this way, the correct modi- 
fied fluctuation spectrum S enters the decay rate. When 
translating this into an appropriate extensive diagram 
(FigEld)), the open-ended single-particle line turns into 
one particle-hole bubble. The latter corresponds to the 
density-density correlator evaluated in the absence of in- 
teractions. In contrast, in a literal (naive) application of 
SP-theory to the many-body problem with screening one 
might be tempted to employ the screened density-density 
propagator (FigfJJe)), which gives incorrect results. 

The source of the difficulty is quite obvious when com- 
paring Fig. [^e) to Fig.^c): The origin of the bubble is 
the single-particle line, which appears since we are inter- 
ested in the coherent propagation of the particle. We are 
not interested in the propagation of a density perturba- 
tion that enters diagram (e). 

The issue discussed here is independent of whether we 
are in the ballistic or the diffusive regime, or of what are 
the details of the model. In any case, the correct applica- 
tion of the SP-theory requires to employ the full potential 
fluctuations (including screening) for S, while looking at 

the bare P}^} , calculated for the non-interacting system. 



VII. SELF CONSISTENT POINT OF VIEW FOR 
THE DEPHASING OF INTERACTING 
ELECTRONS 

We now want to clarify within the SP approach how 
screening should be handled, and demonstrate the consis- 
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tency with the diagrammatic perspective of Section IVIII 
In a dirty metal the motion of electrons is diffusive, and 
there is a screening effect due to the long range Coulomb 
interaction. The fluctuations of the electrostatic poten- 
tial U reflect the many-body fluctuations in the density 
of the electrons via the Coulomb law 



1 

1? 



(37) 



with S and P defined by Eqs. ([7]) and © respectively. 
The equilibrium fluctuations of the induced electrostatic 
fluctuations S" l {q,uj) are trivially related to the dynamic 
form factor P°' , (q,Lu) by the above relation, and can 
be determined self consistently using the fluctuation- 
dissipation relation (see App lCj) : 



1 



vDq 2 



2u) 



1 



(38) 



If we want to calculate the single-particle dephasing rate 
Tp\ then one option is to obtain pM(g,u;) from a blend 
of semiclassical and many-body considerations as in Sec- 
tion IIVI But we are trying in this paper to explore an 
alternative route, where the single particle dephasing rate 
is defined as the many-body dephasing rate Tip per ef- 
fective particle as in Eg . (f3"3"|) . The main difficulty that 
immediately arises, once interacting electrons are con- 
cerned, is how to define the "system" for which the cal- 
culation is done. The proper formulation of the system- 
bath paradigm becomes tricky once electrons are both 
the "system" and the "bath" (see illustration in Figj3j . 
In particular we have to determine what is pl^l (q, to) for 
the case that screening is important. Here we encounter 
a conflict between two opposing points of view, which we 
discuss below. 

There are two ways to determine P^(q,u>). On the 
one hand, within the framework of the heuristic semi- 
classical approach of Section IIVI the spectral function 
PW(q,w) should be given by Eq. (fT5|) which treats the 
electrons as non-interacting. On the other hand, within 
the formal framework it describes the fluctuations of the 
many body density of the electrons, and therefore should 
equal P eq (g,w) as defined after Eq. (|37f . But then one 
realizes that for an interacting system P^(q, oj) is very 
different from P cq (g,tj). So it seems that we have here 
an inconsistency. 

In order to resolve the apparent inconsistency we have 
to clarify the physical meaning of the phrase "dephasing 
rate per particle". This notion is not problematic con- 
ceptually if the environment is a distinct entity (phonon 
bath) . But in the case of a dirty metal this distinction is 
blurred: obviously we cannot regard the same particles 
as both the "system" and "environment" . 

It turns out that a reconciliation of the formal 
fluctuation-dissipation analysis with the heuristic ap- 
proach is possible provided we define the "system" as a 
bunch of non-interacting test particles^, whose density 



fluctuations are simply those of diffusing, non-interacting 
electrons. We can obtain the power-spectrum of the 'mi- 
nority' test particles, starting from Eq. (|C2j) . but replac- 
ing t/totai by the fluctuating potential U produced by the 
'majority' particles. This effectively omits screening ef- 
fects in the calculation of the test particle density fluctu- 
ations: 



(39) 



Note that the p of the system (minority fraction of elec- 
trons) is assumed to be much smaller than the total elec- 
tronic density p clct that appears in the FDT derivation, 
hence one can neglect the back-reaction effect. In other 
words, as far as the minority p is concerned we do not 
take the screening into account, and treat them as non- 
interacting. From Eq.®, together with Eqs.® and © 
and the Einstein relation a = e 2 Dv, we get immediately 



P [N \^) 



L" 



vDq 2 



iu> — Dq 2 



S°«(q,w) (40) 



where for the second line we recalled Eq. JTU) for p(°) 
and Eq. jC8]) for S^. This result agrees with Eq.fH} 
as if we could regard the electrons as diffusing but non- 
interacting. It is important to appreciate that we have 
obtained here a non-trivial profound relation that by- 
passes the heuristic approach which is required in order to 
adopt Eq. (fT5f for strongly interacting electrons: strictly 
speaking the derivation that leads to Eq. (fl"5)) is not appli- 
cable here. Still we get Eq. (j40|) which agrees with Eq. (fT5|) 
by extending the common fluctuation-dissipation reason- 
ing, without the need to introduce a blend of semiclassics 
with Pauli exclusion factors. 



VIII. CONCLUSIONS 

We have presented a straightforward approach to the 
calculation of the dephasing rate within the framework 
of Fermi golden rule picture, and applied it to a many- 
fcrmion system. Starting from the quantum spectra of 
the environment (S) and the system (P), the approach 
(termed here "SP-theory") yields the dephasing rate as 
an integral over the frequency transferred between sys- 
tem and environment during interaction processes. In 
the present paper, we have gone beyond previous at- 
tempts, and considered a full many-fermion system. We 
have argued that this yields results which automatically 
incorporate the crucial physics of Pauli blocking that 
serves to suppress decoherence at low temperatures. The 
many-body dephasing rate can be identified as the single- 
particle dephasing rate times the effective number of 
thermally excited particles susceptible to scattering. The 
use of non-symmetrized spectral functions provides in a 
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natural way the proper cutoff scheme for the frequency 
integral of the " SP formula" , without having to invoke 
the ad hoc cutoff schemes that had been introduced in 
previous publications. 

By defining the single particle dephasing as the many- 
body dephasing per particle one can bypass the need to 
use a blend of semiclassics with Pauli exclusion factors. 
This point of view also provides a natural bridge to the 
diagrammatic approach. Indeed we have shown how the 
results of the SP-theory can be interpreted in terms of di- 
agrams. This has allowed us to address another question, 
namely how SP-theory should be applied in a situation 
in which the system-environment coupling is no longer 
weak. That is the situation relevant for electrons mov- 
ing in a disordered metal, where screening is crucial for 
the structure of the Nyquist noise. We have shown that 
SP-theory should incorporate in such a case the full envi- 
ronment spectrum, alongside the non-interacting density 
spectrum of the system, in agreement with the diagram- 
matic approach. 
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APPENDIX A: THE SEMICLASSICAL PICTURE 
OF DEPHASING 

The notion of dephasing naturally arises in the analysis 
of transport where, loosely speaking, one is interested in 
calculating the probability of a particle to get from one 
point to a different point. Consequently the most popular 
definition of the dephasing factor is based on a semiclas- 
sical picture. Using the Feynman- Vernon formalism, and 
adopting the notations as in Ref.0IH, the propagator of 
the reduced probability matrix in the presence of a ther- 
mal bath, is expressed as a sum over pairs of classical 
trajectories: 



(/ b 



A a A% exp - 



S N [x a ,x b ] 



n 2 



exp 



^[xJ-Sjxb] )(A1) 



The dephasing factor Ptp{t) is defined as a number 
within [0, 1] that characterizes the suppression of the off- 
diagonal elements: One observes that after time t the 
interference contribution, that comes from the off diago- 
nal terms in the double sum, is suppressed by a factor 



where x is the preparation of the bath. In order not 
to complicate the notations, the canonical average over 
the x states is implicit. The unitary operator U[x] gen- 
erates the evolution of the bath given that the particle 
goes along the trajectory x(t). The action Sn[xo,, Xb] is a 
double time integral. Using manipulation as in Refs. 
one obtains the "SP formula" with the symmetrized ver- 
sion of S(q,uj), and the symmetric classical version of 
P(q,uj), namely of Eq.([l4]). 

The semiclassical expression is definitely wrong for 
short range scattering at low temperatures 0HJ , because 
it does not reflect that closed channels cannot be excited. 
This problem with the semiclassical (stationary phase) 
approximation is well known in the theory of inelastic 
scattering. 



APPENDIX B: THE PURITY BASED 
DEFINITION OF DEPHASING 

The following appendix follows the presentation of 
Ref. 12T[T3 |. The natural definition for the dephasing 
factor P v (t) is related to the purity trace(p 2 ) of the re- 
duced probability matrix. Given that the state of the 
system including the environment is ty pn , where p and n 
label the basis states of the particle and the bath respec- 
tively, the full probability matrix is ^pn^p' n i , while the 
purity of the reduce probability matrix is 



= \/trace(p2 y J = 1 /trace( / 3 o 2 n 



\ \f> , ,iLr*„ ,vl/ „ ,,\Si*, ,, 
/ j *p'n' *p"n' *p"n" *p'n" 

p / p"n'n" 



1/2 



.(Bl) 



Assuming a factorized initial preparation as in the con- 
ventional Feynman- Vernon formalism, we propose the 
rate of loss of purity as a measure for decoherence. A 
standard reservation applies: initial transients during 
which the system gets "dressed" by the environment 
should be ignored as these reflect renormalizations due 
to the interactions with the high frequency modes. Other 
choices of initial state might involve different transients, 
while the later slow approach to equilibrium should be in- 
dependent of these transients. In any case the reasoning 
here is not much different from the usual ideology of the 
Fermi golden rule, which is used with similar restrictions 
to calculate transition rates between levels. 

Writing the initial preparation as ^ p °n — Sp, Pa 6 n ,n , 
and using leading order perturbation theory, we can re- 
late P v to the probabilities P t (p,n\po,n ) = "fjml 2 to 
have a transition from the state \po, no) to the state \p, n) 
after time t. The derivation is detailed in Appendix E of 
Ref. [Hill. 

One obtains the result 



P v (t) = P t {po,n Q \p ,n ) 

+ Ptip^Po^olPo^o) 



(B2) 

P t {po,n^n \p ,n ). 



_ e -S N [x A .x u \ = 



v The notation p ^ po or n ^ no implies a summation 

E^Jx C^I^blx ) j (A2) J2p=£ P or Enin,,' respectively. The actual calculation 
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of Pt(p, n\po, no) can be done using Fermi's golden rule 
(FGR) as discussed in the main text. 

Thus we see that within the FGR framework, the pu- 
rity is simply the probability that either the system or 
the bath do not make a transition. Accordingly P v (t) is 
essentially the same as the survival probability P(t) of 
the initial state (the first term in Eq. (IB2[) h In typical 
circumstances the difference between P v (i) and P(t) has 
zero measure weight in the dqduj integration and there- 
fore T v can be identified with the Wigner decay rate of 
system excitations. 



APPENDIX C: THE 
FLUCTUATION-DISSIPATION RELATION 



Note that 
Im 



{Dq 2 + Ana) 2 



us 

47TCT 



(C4) 



The interaction between the electrons and an external 
electrostatic field is described by 7i axt — U ext p elct which 
can be also written as 7i oxt = p axt U . The fluctuation dis- 
sipation relation expresses S oq (<7,a;) using the response 
function a{q,uj) that relates U to — p oxt which is 



Ane 2 



Namely, 



1 



e(q,u) 



(C5) 



The dielectric constant of a metal is defined via the 
linear relation between the total electrostatic potential 
C/ total and an external test charge density p oxt 



1 /47re 2 



P,r- 



(Cl) 



For simplicity we relate here and below to one 
component q of the fields. The total electro- 
static potential is the sum of the external potential 
U cxt = (47re 2 /g 2 )p cxt , and the induced potential U = 
(4ire 2 /q 2 )p Blot , where p clct is the total density of the 
electrons. The dielectric constant can be deduced 
from the equations of motion dp B \ ct /dt = — VJ with 
J = — (cr/e 2 )Vr/ total — £>Vp olGt that leads to the relation 



We 2 )? 2 



relet . j-. 9 

iuj — L)q z 

and hence to U totBl = (l/e)t/ cxt , where 

47rer 



e(q,u) = 1 



iu> — Dq 2 



(C2) 



(C3) 



S*«(q,u>) 



Im 



1 - e~«V T 



leading to 



S*(g,w) 



1 



2w 



o/T 



(C6) 



(C7) 



a q* \ 1 — e 

The Ohmic behavior is cut-off by \u>\ < l/r c and \q\ < 
Xjl where i = v f t c is the elastic mean free path, and 
i> F is the Fermi velocity. Recalling the Einstein relation 
a — e 2 i/D, where v = S^ 1 / L d is the density of states per 
unit volume, we can write this result more conveniently 
as follows: 



1 



2uj 



vDq 2 \\ -e~"/ T 



(C8) 



Note that the electron charge e cancels out from this final 
result for the Nyquist noise spectrum. This well-known 
fact is due to the effects of screening: A larger value of the 
charge would be canceled by a correspondingly stronger 
suppression of density fluctuations. 
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For the purpose of this illustration we assume that the 
velocity and hence the diffusion coefficient of the electron 
do not depend on its kinetic energy (a linear dispersion 
relation). The electrons are added one by one, and at low 
densities form a Boltzmann gas. Only when the chemical 
potential exceeds T one starts to observe the implications 
of Fermi statistics. 

To label one or several electrons as test particles is of course 
meaningless if we take their indistinguishablity into ac- 
count. So the dashed line has some meaning only in a vague 
semiclassical context which we try to bypass. If we literally 
divided the electron into two sub-volumes, the purity loss 
and hence the many-body dephasing rate would be pro- 
portional to the area of dividing surface that separate the 
two spatial regions. 

The notion of "test particles" acquires a more concrete 
meaning within the context of a transport experiment, 
where the "test particles" are the ones contributing to the 
current, or even more to the point, in the context to weak 
localization, to the back-scattered current. 



